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Time  dependent  calculations  of  minor  species  distributions 
require  as  input  the  solar  radiation  field  as  a  function  of  altitude 
and  solar  zenith  angle.  An  isotropic  spherical  model  of  the  radiation  1 

field  has  been  developed  to  determine  the  solar  radiation  field  in  j 

twilight.  Comparison  of  the  spherical  model  with  a  plane  parallel 
model  of  the  multiple  scattering  of  radiation  at  twilight  shows  that:  j 
(1)  for  solar  zenith  angles  less  than  95°,  plane  parallel  solutions 
are  suitable  if  the  initial  deposition  of  solar  energy  is  calculated  1 

for  a  spherical  atmosphere;  (2)  for  solar  zenith  angles  greater  than 

95°,  the  plane  parallel  radiation  field  is  several  orders  of  magnitude 
smaller  than  that  calculated  with  the  spherical  model;  (3)  at 
altitudes  above  40  km,  the  spherical  model  predicts  10  -  20  %  less 
radiation  than  the  radiation  field  calculated  with  the  plane  parallel  J 
model. 

"  \ 

\ 


17.  Key  Words 

Troposphere,  Stratosphere, 
Solar  Radiation,  Twilight, 
Multiple  Scattering 


18.  Distribution  Stotement 

This  document  is  available  to 
public  through  the  National 
Technical  Information  Service 
Springfield,  Virginia  22161 


the 


19.  Security  Clossil.  (ol  this  report) 

UNCLASSIFIED 


20,  Security  Clossil.  (©t  this  page) 

UNCLASSIFIED 


21.  No.  of  Pages 

17 


22.  Price 


"or 


the  solar  radiation  field  at  twilight 


BELOW  100  Kit:  A  SPHERICAL  MODEL 


by 

Donald  F..  Anderson,  Jr. 

E.  0.  Hulburt  Center  for  Space  Research 
Naval  Research  Laboratory 
Washington,  DC  20375 


Abstract 


i 


Time  dependent  calculations  of  minor  species  distributions  require 
as  input  the  solar  radiation  field  as  a  function  of  altitude  and  solar 
zenith  angle.  An  isotropic  spherical  model  of  the  radiation  field  has 
been  developed  to  determine  the  solar  radiation  field  in  twilight. 

Comparison  of  the  spherical  model  with  a  plane  parallel  model  of  the 
multiple  scattering  of  radiation  at  twilight  shows  that:  (1)  for  solar 
zenith  angles  less  than  95°,  plane  parallel  solutions  are  suitable  if  the 
initial  deposition  of  solar  energy  is  calculated  for  a  spherical  atmosphere; 
(2)  for  solar  zenith  angles  greater  than  95°,  the  plane  parallel  radiation 
field  is  several  orders  of  magnitude  smaller  than  that  calculated  with 
the  spherical  model;  (3)  at  altitudes  above  40  kir.,  the  spherical  model 
predicts  10  -  20  %  less  radiation  than  the  radiation  field  calculated 
with  the  plane  parallel  model. 
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I.  Introduction 


Time  dependent  calculations  of  minor  species  distributions  require 
as  input  the  solar  radiation  field  as  a  function  of  altitude  z  and  solar 
zenith  angle  0O.  There  has  been  some  concern  that  results  from  plane 
parallel  models  may  be  incorrect  at  twilight.  In  the  present  study,  an 
isotropic  spherical  model  of  the  radiation  field  has  been  developed  to 
address  the  twilight  problem.  The  model  does  not  include  refraction 
effects,  ground  albedo,  or  aerosols,  but  does  include  molecular  oxygen 
(O2)  and  ozone  (O3)  absorption,  and  multiple  scattering  in  an  inhomogeneous, 
spherical  atmosphere.  Ground  albedo  is  relatively  unimportant  in  twi  light; 
refraction  effects  may  be  neglected  except  when  making  comparison  to 
high  resolution  measurements;  and  the  presence  of  aerosols  slightly 
enhance  the  degree  of  multiple  scattering  at  all  wavelengths.  These 
effects  may  be  included  at  a  later  time  but  the  basic  results  described 
here  will  not  be  affected. 

II.  The  Models 

The  model  atmosphere  employed  is  the  U.S.  Standard  Atmosphere  (1976) 
with  the  O3  profile  from  Nicolet  (1978).  The  O3  cross  section  between 
200  and  800  nm  is  also  from  Nicolet  (1978),  as  is  the  O2  absorption  cross 
section.  The  model  atmosphere  is  shown  in  Figure  1  and  the  O3  and  O2 
cross  sections  in  Figure  2. 

The  spherical  radiative  transfer  model  is  patterned  after  the  model 
developed  by  Anderson  and  Hord  (1977)  for  calculation  of  resonance 
scattering  by  atomic  hydrogen  of  the  solar  Lyman-alpha  flux  incident  on 
planetary  thermospheres  and  exospheres.  The  basic  physical  model  for 
Rayleigh  and  Mie  scattering  is  described  for  the  plane  parallel  case  by 
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Anderson  and  Meier  (1979)  and  Meier  et  al.,  (1982).  Detailed  nunerical 
results  for  the  wavelength  region  200-800  nm  in  the  troposphere  and 
stratosphere  are  described  by  Nicolet  et  al .  ,  (1982). 

Briefly,  both  the  plane  parallel  and  spherical  theory  used  in  this 
study  solve  the  integral  equation  of  radiative  transfer  for  the  flux  into  a 
volume  element  F.  F  is  normalized  to  the  solar  flux  at  the  top  of  the 
atmosphere  F^ ,  where  F^  is  in  units  photons  cm^s-ln^-l.  The 
plane  parallel  model  assumes  azimuthal  symmetry  and  utilizes  Chandrasekhar's 
(1960)  E  functions  to  reduce  the  problem  to  one  dimension:  namely  F 
as  function  of  vertical  optical  depth  t  ,  or  altitude.  In  the  spherical 
code,  azimuthal  symmetry  is  assumed,  but  F  is  now  an  interlocking  function 
of  t  and  0O.  Thus,  on  multiple  scattering,  photons  are  allowed  to 
migrate  in  both  t  and  0O.  It  will  be  shown  below  that  the  effects 
of  this  coupling  can  be  very  large  in  twilight. 

III.  Results 
A.  Calculation  of  FQ 

When  plane  parallel  models  are  employed  to  calculate  F(  t  ,  0o),  the 
attenuated  direct  flux  F0(  t  ,  0o)  is  usually  given  by 

Fo(  T  »  0O)  =  exPf-(  T  +  Tp)  sec  0O]  (1) 

t  is  related  in  a  straightforward  way  to  z  noting  that 

T  “  °s  /  N(z)<te  (2) 

where  og  is  the  cross  section  for  Rayleigh  scattering,  and  N  is  the  neutral 
atmosphere  density  at  altitude  z  in  the  atmosphere.  An  absorption  which 
removes  a  photon  from  further  scattering  is  designated  a  pure  absorption, 
and  the  pure  absorption  optical  depth  is 

tp  -  /  l  0O3nO3(z>  +  ao2No2(z))dz 


(3) 


where  the  o's  are  the  pure  absorption  cross  sections.  If  the  pure 
absorber  is  homogeneously  mixed  with  the  scatters,  then  Chapman  functions 
(Smith  and  Smith,  1972)  may  be  used  in  place  of  sec  0Q  in  (1)  to  account  for 
effects  on  F0  of  sphericity.  Since  O3  is  not  homogeneously  mixed  neither 
Chapman  functions,  nor  sec  0o  is  appropriate.  Instead,  the  true  slant 
optical  depth  must  be  calculated  taking  into  account  the  inhomogenit ies  due 
to  ozone.  A  comparison  of  (1)  with  the  correct  F0  is  shown  in  Figure  3 
at  0Q  =  80,  86,  and  88°,  for  A  =  325  nm ,  ip  =  0.14  and  A  =  325  nm, 

Tp  =  0,  respectively.  Clearly,  below  50  km,  spherical  effects  are 
important  for  0Q  >  80°.  No  comparisons  are  made  for  ©Q  90, 
since  the  plane  parallel  F0  is  not  defined.  In  Figure  4,  F0  is  shown 
for  A  =  325  nm  and  0Q  from  0°  to  98°.  The  wavelength  A  =  325  nm 
is  chosen  throughout  for  illustrative  purposes,  since  both  pure  absorption 
and  multiple  scattering  effects  are  important.  However,  the  model  does 
cover  the  entire  range  200  -  800  nm. 

B.  Calculation  of  F 

Values  of  F  at  A  =  325  nm  are  shown  in  Figure  5  for  0Q  =  60°  to 
104°.  The  dashed  curves  are  the  values  of  F  calculated  with  the  correct 
spherical  F0 ,  but  the  multiple  scattering  is  assumed  to  take  place  in  a 
plane  parallel  atmosphere.  It  is  clear  from  Figure  5  that  for  0O  <  92° 
the  plane  parallel  solution  with  a  spherical  F0  is  adequate  to  describe 
the  multiple  scattering.  However,  for  ©0  >  96°  the  plane  parallel 
solutions  are  orders  of  magnitude  lower  than  the  spherical  solutions  and 
are  off  scale. 

Another  result,  not  apparent  on  the  log  plot  in  Figure  5,  is  that 
above  40  km,  F  (spherical)  falls  about  10  to  20%  below  F  (plane  parallel) 
at  all  0O.  This  result  has  been  noted  before  (Strickland  and  Anderson, 


1973)  and  occurs  because  a  photon  emitted  at  a  given  altitude  in  a  spherical 
atmosphere  has  a  greater  probability  of  escape  from  the  medium  than  in  a 
plane  parallel  atmosphere.  As  z  ->  “  the  escape  probability  approaches 
0.5  in  a  plane  parallel  medium,  and  1.0  in  a  spherical  medium. 

As  an  example  of  the  wavelength  dependence.  Figures  6  and  7  show 
F(0O)  for  X  =  250  and  450  ran,  respectively.  At  short  wavelengths 
(  X  <  325  nm)  ozone  absorbs  much  of  the  radiation  so  that  F  is  low,  and 
is  essentially  zero  below  the  ozone  maximum.  As  X  increases  beyond  350  nm 
F  gradually  decreases  in  proportion  to  the  Rayleigh  optical  depth,  and 
thus  the  multiple  scattering  decreases. 

IV.  Discussion 

The  purpose  of  this  research  effort  was:  (1)  to  develop  a  spherical 
model  of  the  multiple  scattering  of  solar  radiation  in  the  twilight 
atmosphere;  and  (2)  to  establish  a  range  of  validity  for  the  use  of  plane  parallel 
models  in  twilight.  The  Initial  model  has  been  developed.  A  more  sophisticated 
model  may  be  developed  in  the  future  which  includes  such  effects  as 
refraction  of  sunlight,  anistropic  scattering  by  aerosols,  ground  albedo, 
and  clouds.  But  the  result  of  comparison  with  the  present  model  and  a 
plane  parallel  model  suggests  that  for  0Q  <  95°,  the  plane  parallel 
solution  is  sufficient,  if  the  proper  F0  is  used.  This  conclusion  is  stated 
with  the  caveat  that  the  spherical  solutions  are  10-20%  lower  than 
the  plane  parallel  solutions  at  all  0o  for  z  >  40  km. 

For  0Q  >  95°  plane  parallel  solutions  are  not  satisfactory,  and 
the  future  improvement  of  the  present  model  depends  to  some  extent  on 
whether  or  not  the  enormous  differences  in  F  between  the  two  models  have 
any  effect  on  the  photochemistry  in  the  atmosphere  below  100  km. 
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Figure  1.  Neutral  Atmosphere  Density  Distribution.  Ozone  is  scaled 
by  a  factor  of  100. 
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Figure  2.  O2  and  O3  pure  absorption  cross  sections.  O2  cross  section  is 
scaled  upward  by  a  factor  of  10. 


Figure  3.  Percent  error  in  plane  parallel  Fn  when  compared  to  spherical 
F0  at  X  =  325  nir.  and  0o  as  indicated,  hashed  and  solid 
curves  are  for  in  =  0  and  0.14  respectively. 


Figure  4.  Spherical  FQ  at  A  «=  325  run  and,  0O  as  indicated 
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